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ABSTRACT: We show that a bounded, isolated quantum system of many particles in 
a specific initial state will approach thermal equilibrium if the energy eigenfunctions which 
are superposed to form that state obey Berry's conjecture. Berry's conjecture is expected 
to hold only if the corresponding classical system is chaotic, and essentially states that the 
energy eigenfunctions behave as if they were gaussian random variables. We review the 
existing evidence, and show that previously neglected effects substantially strengthen the 
case for Berry's conjecture. We study a rarefied hard-sphere gas as an explicit example of 
a many-body system which is known to be classically chaotic, and show that an energy 
eigenstate which obeys Berry's conjecture predicts a Maxwell-Boltzmann, Bose-Einstein, 
or Fermi-Dirac distribution for the momentum of each constituent particle, depending on 
whether the wave functions are taken to be nonsymmetric, completely symmetric, or com- 
pletely antisymmetric functions of the positions of the particles. We call this phenomenon 
eigenstate thermalization. We show that a generic initial state will approach thermal equi- 
librium at least as fast as 0(%/ A)t _1 , where A is the uncertainty in the total energy of 
the gas. This result holds for an individual initial state; in contrast to the classical theory, 
no averaging over an ensemble of initial states is needed. We argue that these results 
constitute a new foundation for quantum statistical mechanics. 
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I. INTRODUCTION 



Take some helium atoms, put them in one corner of a well insulated box, and let them 
go. Wait a while, then punch a small hole in the side of the box. As the atoms emerge, one 
by one, measure their momenta. Make a histogram, plotting the fraction of atoms with 
magnitude of momentum between p and p + dp. 

Every physicist knows what the result of this experiment will be. The histogram will 
be very well approximated by the Maxwell-Boltzmann distribution 

/ MB (p,T) = (27vmkT)- 3 / 2 e-^ 2mkT (1.1) 

multiplied by 4nFp 2 dp, where F = (ir/SmkT^^p is a flux factor. Here m is the mass of 
a helium atom, k is Boltzmann's constant, and T is the "temperature," a number which 
will depend on how the atoms were originally put into the corner of the box, how much 
space they occupied, and other details of the initial conditions. The challenge is to derive 
this result from first principles. 

The biggest problem with a theoretical analysis of this particular experiment is the 
need to treat the hole in the box in a reasonable way. It is much easier to study the 
following thought experiment instead. Suppose, after preparing the system in its initial 
state, we are able to measure the momentum of one atom at a specific (but arbitrary) time 
t. Suppose further that, after having made this measurement, we can empty out the box, 
and then start it off again with the system in exactly the same initial state. We now do 
this repeatedly, each time measuring the momentum of one atom after exactly the same 
amount of time t has passed. We make a histogram of the results. 

Let us analyze this experiment, beginning with classical mechanics as the underlying 
theory. We take the hamiltonian for N atoms to be 

i=l i<j 

where we take V(r) to be a hard-sphere potential: 

f +00 for r < 2a 
V(r) = (1.3) 
I for r > 2a . 

We assume perfectly reflecting boundary conditions at the walls of the box. The atoms 
initially have some definite total energy U. The phase space of this system is known to be 
fully chaotic, with no invariant tori for any value of U [1] . Thus the motion in phase space 
on any constant energy surface is ergodic and mixing. (For a review of classical chaos 
theory, see, e.g., [2-4].) 
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However, this is entirely irrelevant if we always start out with exactly the same initial 
state, and always make the measurement after exactly the same amount of time has elapsed. 
The momentum of the measured atom (assuming that it is always the same atom) is 
determined exactly by the initial conditions, and so will always be the same. To have 
any hope of getting a distribution of momenta, we must average over either the initial 
conditions or the times of measurement or both. 

If we keep the initial conditions fixed, ergodicity implies that the system wanders all 
over the constant energy surface. (This assumes that we have not started the system off at 
a point located on a periodic orbit; such points form a set of measure zero.) If we divide 
the constant energy surface into many patches of equal area (and that area is not too 
small) then after a certain finite time the system will, to a very good approximation, be 
equally likely to be in any one of these patches at any later time. Conversely, if we permit 
a range of initial conditions, mixing implies that, if the measurement time is fixed but 
not too early, the system will once again, to a very good approximation, be equally likely 
to be in any one of the equal-area patches. The rule that equal-area patches are equally 
likely is just the usual formulation of the microcanonical ensemble, which after a little 
work leads to eq. (1.1) for the fraction of atoms with momentum in a range of d 3 p around 
p, with the temperature T simply given by the ideal gas formula U = ^NkT. In short, 
if we do a modest average over either the initial conditions or the times of measurement, 
then classical chaos results in classical thermalization. (For an elementary review of chaos 
theory as it applies to classical statistical mechanics, see [5].) 

On the other hand, if we have a weakly perturbed integrable system (for example, 
harmonic oscillators with small nonlinear couplings), then according to the Kolmogorov- 
Arnol'd-Moser theorem [6], its phase space is foliated by invariant tori almost everywhere, 
and we do not expect it to thermalize. If it is partially integrable, with some invariant tori 
embedded in an otherwise chaotic phase space, then the system may or may not thermalize, 
depending on the initial state. 

These results from classical mechanics are clear and powerful, and provide a satisfying 
explanation of statistical behavior in classical systems which exhibit chaos. However, 
we know that the real world is ultimately described by quantum mechanics, and we so 
should seek the quantum analog of the classical analysis. We would like to know, for 
example, what property a quantum system must possess, analogous to classical chaos, so 
that "most" of its initial states thermalize, in the sense discussed above. Furthermore if 
a quantum system does possess this property (whatever it may be), then we might hope 
that the inherent uncertainties in quantum mechanics lead to a thermal distribution for the 
momentum of a single atom, even if we always start with exactly the same initial state, and 
make the measurement at exactly the same time. If this is true, then quantum mechanics 
automatically provides the "coarse graining" which is missing [3,7] in the classical theory. 

I will argue that the property needed for thermalization of a quantum system is the 
validity of Berry's conjecture [8-10]. For a quantum gas of hard spheres, Berry's conjecture 
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states that each energy eigenfunction appears to be a superposition of plane waves (in the 
3iV dimensional coordinate space) with random phase and gaussian random amplitude, 
but fixed wavelength. In general, Berry's conjecture is expected to hold only for systems 
which exhibit classical chaos in all or at least most of the classical phase space. As already 
noted, a hard sphere gas meets this condition. 

We will see that Berry's conjecture leads to either Bose-Einstein, Fermi-Dirac, or 
Maxwell-Boltzmann statistics, depending on whether the wave functions are chosen to 
be completely symmetric, completely antisymmetric, or nonsymmetric functions of the 
positions of the iV atoms. Furthermore we will find that any nonthermal features of the 
initial distribution of momenta decay away at least as fast as 0(/i/A)t _1 , where h is 
Planck's constant and A is the uncertainty in the total energy. Thermal behavior thus 
appears for a very wide range of possible initial states, without assuming that the system 
interacts with an external heat bath, or any other environmental variables. We also do 
not need to take any averages over initial states, times of measurement, or hamiltonians, 
or make any unjustifiable approximations to the quantum equation of motion, such as 
truncation of the BBGKY hierarchy. 

The rest of this paper is organized as follows. In Sect. II, we review Berry's conjec- 
ture for a system of hard spheres in a box. In Sect. Ill, we show how Berry's conjecture 
leads, in the limit of low density and high energy, to a Maxwell-Boltzmann distribution 
for the momentum of a single atom in the gas, with a temperature that is related to the 
total energy by the ideal gas law. In this section we treat the atoms as distinguishable, 
making no assumptions about the symmetry of the wave function of the gas under ex- 
change of individual atoms. In Sect. IV, we evaluate the effects of certain corrections 
to Berry's conjecture known as "scars," and reconsider some of the numerical results on 
Berry's conjecture which have appeared in the literature. In Sect. V, we examine time evo- 
lution beginning with a nonthermal initial state, and study the approach to equilibrium. 
In Sect. VI, we consider wave functions which are completely symmetric or antisymmetric 
under exchange of atoms, and see that these lead to Bose-Einstein or Fermi-Dirac dis- 
tributions for individual momenta, respectively. Conclusions, speculations, and possible 
extensions are presented in Sect. VII. 



Consider a system of N hard spheres, each of radius a, in a cubic box with edge length 
L + 2a. Call the energy eigenvalues U a and the corresponding eigenfunctions ip a (X), where 
X = (x l5 . . . , xtv) denotes the 3iV coordinates, and P = (pi, . . . , pjv) will denote the 3iV 
conjugate momenta. We take the wave functions to be defined on the domain 



II. BERRY'S CONJECTURE 




Xi,...,Xjv 



\L < x a ,2,3 < +\L; |xj - x 



> 2a | , 



(2.1) 
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with the boundary condition that each ^ a (X) vanish on the boundary of D. For now we 
assume that ip a (K) has no symmetries under exchange of individual Xj. 

The energy eigenfunctions tp a (X.) can always be chosen to be everywhere real, and 
can be written as 

^«(X) = M a J d 3N P A a (P) 5(P 2 - 2mU a ) exp(zP-X/ft) , (2.2) 
where J\f a is a constant to be determined by the normalization condition 

/ «TV£(X) = 1, (2.3) 

JD 

and where A* a (P) = A a (— P). For this system, Berry's conjecture is equivalent to assuming 
that A a (P) can be treated as a gaussian random variable with a two-point correlation 
function given by 

<A a (P)A /3 (P')) EE = 5 aP 5 3N (P + P')/5(P 2 - P' 2 ) . (2.4) 

Here 5 3Ar (P) is the 3iV-dimensional Dirac delta function and S(x) is the one-dimensional 
Dirac delta function. The subscript EE stands for "eigenstate ensemble." This is a ficti- 
tious ensemble which describes the properties of a typical energy eigenfunction. Individual 
eigenfunctions behave as if they were selected at random from the eigenstate ensemble. 
Berry's conjecture also asserts that the eigenstate ensemble is gaussian, so that all multi- 
point correlation functions are given in terms of the two-point correlation function; e.g., 

<A a (P 1 )^(P 2 )A 7 (P 3 )^(P 4 )) EE = <A,(Pi)MP 2 )> EE <A 7 (P 3 )A 5 (P 4 )) ee 

+ (,4 Q (Pi)A(P3)> EE <^(P 4 )^(P 2 )> EE 
+ <A a (P 1 )A 5 (P 4 )) EE <^(P 2 )A 7 (P 3 )) EE . (2.5) 

Of course, each A a (P) must give back a ip a (X-) which vanishes on the boundary of D; this is 
not a stringent requirement on A a (P) at high energy, where ifj a (X) has many wavelengths 
between any two segments of the boundary. We will say more about the requirement 
of high energy shortly. Meanwhile, for a more general but less transparent definition of 
Berry's conjecture, see sect. 4; for related mathematical results, see [11,12]. 

Berry's conjecture is based on semiclassical reasoning, and is manifestly untrue for 
systems whose classical phase space is foliated almost everywhere by invariant tori [13]. It 
has been investigated numerically for simple systems which are fully chaotic classically, such 
as a single particle in a two-dimensional stadium-shaped box [14-18], or a single particle 
on a two-dimensional surface with constant negative curvature and periodic boundary 
conditions [19,20]. In these systems, Berry's conjecture is found to be valid for eigenstates 
of sufficiently high energy, and its validity has even been suggested as a good definition of 
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chaos in a quantum system [21]. However, an important caveat is the existence of "scars" 
on some energy eigenfunctions, regions of enhanced value of if) a (X) which follow the paths 
of the most stable classical periodic orbits [14,15]. For now we will ignore the scars, since 
their presence will not alter any of our conclusions. We will discuss them in more detail in 
sect. 4. Also, there we will argue that some numerical results which have been interpreted 
as evidence against Berry's conjecture actually provide evidence for it. 

Even in a system which is fully chaotic classically (like a hard-sphere gas), Berry's 
conjecture will certainly not be valid for eigenfunctions which are too low in energy. The 
low-lying states necessarily have specific structure: the ground state, for example, is node- 
less. A rough criterion for the validity of Berry's conjecture is that the average wavelength 
of each particle be small enough to "see" the features which produce classical chaos [16]. 
For the hard-sphere gas, the relevant feature is the nonzero radius a of each particle. Antic- 
ipating a bit and defining a temperature T a for each energy eigenvalue U a via the ideal gas 
formula U a — ^NkT a , and further defining a "thermal" wavelength A Q = (2nfi 2 /mkT Q ,) 1 / 2 , 
then the criterion for the validity of Berry's conjecture is A a < a. Numerically, this be- 
comes T a > (300 /ma 2 ) Kelvin, where a is in angstroms and m is in amu. 

It turns out that getting explicit results will also require us to work at low density, 
Na 3 <C L 3 . Combining this with A Q < a, we see that we need to have A^A 3 C I 3 , a 
condition which is also required, in quantum statistical mechanics, for the Bose-Einstein 
and Fermi-Dirac distributions to be well approximated by the Maxwell-Boltzmann distri- 
bution. 

Let us now consider the eigenfunctions in momentum space: 

^ Q (P) = h~ 3N / 2 [ d 3N X ^a(X) exp(-zP-X/ft) 
Jd 

/+oo 
d 3N K A a (K)5(K 2 - 2mU a ) 8 3 D N {K - P) , (2.6) 
-oo 

where we have defined 

S 3 D N (K) = h~ 3N [ d 3N X exp(zK-X/ft) . (2.7) 
Jd 

If the condition needed for Berry's conjecture is satisfied (A a < a), and we are in the low 
density regime (Na 3 <^ L 3 ), then we can make the substitutions 

6% N (0) - (L/h) 3N , 

5 3 D N (P) - 5 3N (P) , (2.8) 
[5 3 D N (V)] 2 ^(L/h) 3N 5 3N (V) . 
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Now using eqs. (2.2), (2.4), (2.6), and (2.8), we find that (for A a < a and Na 3 < L 3 ), 

<C(P)^(P')> EE = 5 aP Ml h 3N 5(V 2 - 2mU a ) 8l N (P - P') , (2.9) 
which will play a key role in the next section. 

III. EIGENSTATE THERM ALIZATION 

Let us now put our gas of N hard spheres into some initial state specified by the 
momentum-space wave function 

^(P,0) = ^C Q ^ a (P). (3.1) 

a 

We take the energy eigenfunctions to be ortho normal, and also assume ip(P, 0) to be 
normalized, so that J2 a |C«| 2 = 1- The expectation value of the energy is then 

^ = ^|Ca| 2 ^a, (3.2) 

a 

and the uncertainty in the energy is A, where 

A 2 = ]T|C a | 2 (t/ Q -£/) 2 . (3.3) 

a 

We will assume that A <^ U. The initial wave function will evolve in time according to 
the Schrodinger equation: 

$(P,t) = ^C Q expHtU/ft)^ Q (P) . (3.4) 

a 

Now return to the thought experiment in which the system is repeatedly prepared in the 
same initial state (specified by the C a 's), and the momentum of one atom is repeatedly 
measured after the same elapsed time t. The theoretical prediction for the fraction of 
atoms with momentum in a range d s p around p is / QM (p, t)d 3 p, where 

/qm(pi,*)= [ d 3 p 2 ...d 3 p N \^(P,t)\ 2 




^QC^e^-W $ a/3 ( pi ) . (3.5) 

a/3 
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In the last line we have introduced 



d> Q/3 ( Pl ) = J d 3 P2 ...d s PN r a (p)i>p(p) , (3.6) 

which obeys the normalization condition 

J dV$ a/3 ( P i) =5 aP . (3.7) 

If the system thermalizes, then after some time has passed, / QM (pi^) should be equal to 
the Maxwell-Boltzmann distribution of eq. (1.1), although (as in the classical case) some 
modest averaging over either the initial conditions (the C a 's) or the times of measurement 
(the value of t) might first be necessary. Furthermore the temperature T should be given 
at least approximately by the ideal gas law U = ^NkT, with a fractional uncertainty of 
order A/U. 

To practice on a simple example, let us study the case where the initial state is a 
single energy eigenstate. This is, of course, unphysical: we cannot actually prepare such 
a state in a time less than 0(K/5), where 5 is the mean energy level spacing near U a [22]. 
This is fantastically small in any realistic case [8 a = l/n Q , where n a is given by eq. (4.6)], 
so that %/5 a is much longer than the age of the universe. Nevertheless taking the initial 
state to be an energy eigenstate will turn out to be an instructive exercise. 

In this case, eq. (3.5) becomes simply f QM (pi,t) = $ QQ (pi), which is independent of 
time. We now study the properties of $ QQ (pi) in the eigenstate ensemble introduced in 
the previous section. Assuming high energy (A Q < a) and low density (Na 3 <C L 3 ), it 
follows from eqs. (2.8), (2.9), and (3.6) that 

($ aa ( Pl )) EE = A/; 2 L 3N J d 3 p 2 . . . d 3 p N S(P 2 - 2mU a ) . (3.8) 



We introduce the useful formula 



I D (x) ee J d D P5(P 2 -x) 

(nx) D / 2 



V(D/2)x ' 

and use it to fix N~ 2 = L 3N I 3N (2mU a ) via eq. (3.7). We then find 

/ih f n \\ hN-z(2mU a -pi) 
<$aa(Pl)> - " h N (2mU a ) 

= T(3N/2) ( 1 \ 3/2 / _ p 2 x( 3jV - 5 )/ 2 
~ r((3AT-3)/2) \2iimUa) \ 2mU a 



(3.9) 



(3.10) 
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If we now set U a = ^NkT a and take the large N limit, we get 

<d> aa ( Pl )) EE = {2-KmkT a )-V 2 e-rt' 2mkT « , (3.11) 

which is precisely f MB (pi,T a ). 

Let us note first that, given eq. (3.8) as a starting point, eqs. (3.9-11) simply recapit- 
ulate a standard derivation of the canonical ensemble from the microcanonical [23] . 

More importantly, we must study the fluctuations of <fr aa (pi) about its average value 
in the eigenstate ensemble. We begin by defining 

[A$ a/3 (Pi)] 2 ee <|^(Pi)| 2 ) EE " |(^(Pi)> E J 2 • (3-12) 
Using eqs. (2.5), (2.8), (2.9), and (3.6), we find 

[A<MPi)] 2 = A/^A/J (Lh) 3N J d 3 p 2 . . . d 3 p N d 3 p' 2 . . . d 3 p' N 

x 5(P 2 - 2mU a ) 5(P' 2 - 2mU p ) 6%*(P - P') . (3.13) 

Before evaluating eq. (3.13) explicitly, we can see that it will be very small: if we replace 
5™(P — P') by its maximum value (L/h) SN everywhere, the right-hand side of eq. (3.13) 
becomes ($ QQ (pi)) (^^(pi)) [cf. eq. (3.8)]. Of course this replacement results in 
a huge overestimate of [A<& a( g( P i)] 2 , since in fact 5™(P — P') is close to zero almost 
everywhere. Thus we will have, in particular, A$ aa (pi) <^ ^$ aa (pi)) EB . Furthermore 

we see why A$ QQ (pi) is so small: if) a (P) has fluctuations of order one in the eigenstate 
ensemble, but these are washed out when we integrate over most of the momenta. (In 
Sect. IV we will see that the same fate befalls the "scars" mentioned in Sect. II.) 

We now turn to the evaluation of eq. (3.13), to find out just how small it is. We will 
need to know a bit more detail about Sj^ (P — P') than the substitution rules of eq. (2.8). 
We therefore approximate it with a gaussian: 

6 3 D N (P - P') ~ (L/h) 3N exp[-(P - P') 2 L 2 /4Trh 2 ] . (3.14) 

In the low density regime (Na 3 L 3 ), using eq. (3.14) instead of eq. (2.7) changes the 
result only by an overall constant of order one. Substituting eq. (3.14) into eq. (3.13), 
setting a = (3, and taking the large iV limit yields 

A* aa ( Pl ) = 0(l)JV^ e -3"/4^ . (3 . 15) 

Since we have L ^> A Q , we see that the fluctuations in $ QQ (pi) about f MB (pi,T a ) are 
negligibly small for large N. That is to say, an energy eigenstate which satisfies Berry's 
conjecture predicts a thermal distribution for the momentum of a single constituent par- 
ticle. We will refer to this remarkable phenomenon as eigenstate thermalization. 
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Given a system which exhibits eigenstate thermalization, it is not hard to understand 
why almost any initial state will thermalize. In fact, the problem now is primarily to 
prevent the system from having a thermal distribution for the momentum of each particle 
at all times. To do so at t = 0, we must carefully superpose energy eigenstates in order 
to construct an initial state with whatever nonthermal features we might want. Once 
this superposition is prepared, however, the delicate phase relationships we have set up 
to avoid thermal behavior will gradually be destroyed by hamiltonian time evolution, and 
the system will equilibrate. We will see how this works in more detail in Sect. V. First, 
however, we digress briefly to discuss the scars, and check to see that they do not change 
any of our conclusions so far. 

IV. FADED SCARS 

The theory of scar formation has been developed by Heller [15,18], Bogomolny [24], 
and Berry [25,10]. We will be rather schematic here; readers unfamiliar with scar theory 
should consult the cited references for more details. 

We first consider any system governed by a hamiltonian H(P, X) which is time-reversal 
invariant and which results in classical chaos. For consistency of notation with the previous 
sections, we take P and X to be vectors with 3iV components. 

We begin by introducing the Wigner density for an eigenstate, 

p a (P,X) = h ~ 3N J d3Ns exp(zP-S/ft) i(; a (X + ^S)i; a (X - iS) , (4.1) 
where ifj a (X.) is real. The Wigner density has the useful properties that 

J d w Pp a (P,X)=^(X), 

(4.2) 

J «Tp Q (P,X) = |^(P)| 2 , 
which imply the normalization condition 

J d 3N P d 3N X p a (P, X) = 1 . (4.3) 
Scar theory begins with a semiclassical formula for p a (P, X): 

p a (P, X) = n" 1 h~ 3N 6{H(P, X) - U a ) {l + A v etSp/h ex P [(*A)Z P - W P (T P ) -Z p ] } . 

v 

(4.4) 
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The sum is over all periodic orbits on the surface with constant energy U a ; S p is the action 
of the orbit; T p is a coordinate in phase space along the orbit; Z p are the 6N — 2 coordinates 
in the energy surface which are perpendicular to the orbit; and A p and W p (T p ) are purely 
classical quantities which depend on the monodromy matrix of the orbit. The constant n a 
is fixed by the normalization condition, eq. (4.3), and can be interpreted physically as the 
energy eigenvalue density near U a . If we ignore the sum over periodic orbits in eq. (4.5), 
we obtain the "Weyl rule" for n a : 



For bosons or fermions, the right-hand side should be divided by Nl. Even so, n a is 
fantastically large in any realistic case [22]. 

The key point for scar theory is that the periodic-orbit terms in eq. (4.4) have no %- 
dependent prefactors; the peak height of each term is controlled by the classical quantity 
A p . A short periodic orbit can have an A p which is greater than one; this produces an 
obvious "scar" in phase space along the path of the orbit. 

We are interested, however, in $ QQ (pi), and so we must integrate p Q (P,X) over all 
3N components of X, and all but three of the 3iV components of P. If we consider isolated 
periodic orbits, we see from eq. (4.4) that an integral over one of the 6N — 2 components 

1 /2 

of Zp yields a prefactor of fi 1 . Thus the contribution of each isolated periodic orbit to 
$ QQ (pi) is suppressed, relative to the leading term, by h^ 6N ~ 3 ^ 2 , which means that each 
individual scar on $ QQ (pi) is totally negligible. 

Nonisolated periodic orbits are a little more complicated, since moving off a non- 
isolated orbit in some directions in phase space merely puts the system onto another 
nonisolated orbit in the same "family" ; there are a finite number of these families. For the 
hard-sphere gas, the nonisolated orbits consist of motions where the spheres bounce off 
the walls but never collide with each other [26] . A given nonisolated orbit of this type can 
in general be deformed into another one by moving in any direction in coordinate space. 
Thus when integrating p Q (P, X) over the 3N components of X, we do not get any factors 
of h 1 / 2 . We do, however, get a net factor of ft( 3iV_3 )/ 2 from integrating over 3N — 3 com- 
ponents of P. Thus, while the contribution of a family of nonisolated orbits to $ QQ (pi) 
is much larger than the contribution of a single isolated orbit, it is still negligible. These 
conclusions are supported by the much more detailed calculation of Gaspard [26] for the 
periodic-orbit corrections to n a for the hard-sphere gas. 

Therefore, in computing $ QQ (pi), we can safely ignore the short isolated orbits and all 
of the nonisolated orbits. The most modern version of Berry's conjecture [10] then assigns 




(4.5) 



In the case of a hard-sphere gas of 3N distinguishable particles, this becomes 




(4.6) 
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the origin of the gaussian fluctuations in ^ a (X) to the long isolated orbits. Thus V'a(X) is 
conjectured to behave like a gaussian random variable with a two-point correlation function 
embodied by the elegant formula [8,11] 

<p a (P, X)) ee = n- 1 h~ 3N 5{H(P, X) - U a ) . (4.7) 

There has been important progress recently [27] in bringing the long isolated orbits under 
analytic control, but so far Berry's conjecture remains just that. Even if a rigorous proof 
is eventually found, it is likely to apply only for asymptotically high energies. We turn, 
therefore, to a discussion of the existing numerical evidence. 

Berry's conjecture has been studied numerically in some two-dimensional systems, 
such as a particle in a stadium-shaped box [14-18]. One popular object to study is the 
correlation function [8] 



C a (s) = J d 2 x ^ a (x+ |s)V>a(x- \s) 



d pd x exp(— ip-s/h) p a (p, x) , (4.8) 

where in the simplest case the integral over x covers the entire box, whose area we will 
call L 2 . For a particular eigenstate, the numerically computed C a (s) is compared to its 
expectation value in the eigenstate ensemble: 

(C«(s)) EE = J d 2 pd 2 x exp(-zp-s/ft) (p a (p,x)) BB 
= n" 1 h~ 2 L 2 J d 2 p 5(p 2 /2m - U a ) 

= Jo(k a s) , (4.9) 

where Jq(z) is a Bessel function, k a = (2mU ce /h 2 ) 1 ^ 2 , and s = |s|. In [16], only moderately 
good agreement was found with this prediction, with discrepancies of approximately 0.1 for 
k a = 65 and L = ^fix; see fig. 7 of [16]. However, these discrepancies are entirely explained 
by consideration of the fluctuations in C a (s) which are predicted by the eigenstate ensemble: 



) 2 

EE 



[AC a (s)] 2 = <^(s)) EE -<C a (s)> 

= J d 2 xd 2 y <^ a (x+ §s)V>a(y + | S )) EE <^«( X - ^s)^ a (y- |s)) 

+ <^ a (x+ I S )^a(y- |s)) ee <^ q (x+ I S )^a(y- I*)) 

= L~ 4 J d 2 x d 2 y Jq (k a \x - y|) + J (/c a |x - y + s|) Jo(/c a |x - y - s|) 



(4.10) 
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The first term in the last line dominates over the second for all s, and is 0(l/k a L). That 
is, we expect discrepancies of roughly (k a V)~ x l 2 between C a (s) as computed numerically 
for a particular eigenstate and (C a (s)) EE as given by eq. (4.9). This is exactly what is 
seen in fig. 7 of [16]. Similar comments apply to figs. 14-17 of [20]. The fact that these 
discrepancies are predicted by the eigenstate ensemble does not seem to have been noticed 
previously. 

We note finally that the gaussian nature of the eigenstate ensemble [which is used 
crucially in eq. (4.10)] has also been directly tested. The prediction is that 

<^ n (x)) EE = (2n - 1)!! <^(x)>; E = (2n - 1)!! L~ 2n , (4.11) 
but this just means that the probability distribution for the amplitude ifj at any point is 

P(V0 = ^exp[-±LV] , (4-12) 
which is well supported by the numerical results [16,19,20]. 



V. TIME EVOLUTION AND EQUILIBRATION 



In Sect. Ill, we saw that an individual energy eigenstate which satisfies Berry's conjec- 
ture predicts a thermal distribution for the momentum of each constituent particle. Now 
we must see what happens when we consider more general states. We will once again 
express the initial state as a wave function in momentum space, and expand it in energy 
eigenfunctions, as in eq. (3.1). We assume that the uncertainty A in the total energy, 
eq. (3.3), is much smaller than the average energy U, eq. (3.2). This is easy to arrange in 
practice. 

The predicted momentum distribution of a single particle at time t is f QM (Pi,t), as 
given by eq. (3.5). Now take the average of eq. (3.5) in the eigenstate ensemble. From 
eq. (2.4), it is immediately obvious that ($ a /3(pi)) = if a ^ f3, and so we get 

(/qm(Pi,*)) bb = l^l 2 < $ ««(Pi)>ee 

a 

= \ C <*\ 2 (2nmkT a )- 3 / 2 e~^/ 2mkT - , (5.1) 

a 

where we have used eq. (3.11) in the second line. Now, since A <C U, we can with negligible 
error replace each T a in this sum with T, where U = |iV/cT, and then \C a \ 2 = 1 gives 
us 

</ Q M(pi,t)> EE = /MB(pi,T)[l + 0(A/t/)] , (5.2) 

the desired result. 
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Once again, though, we must study the fluctuations of / QM (pi>*) that are predicted 
by the eigenstate ensemble. We can write 

/qm(pi,*) =/mb( P i,T) +J2 C «Cp^ (Ua - U " )t/h ®MPi) , (5-3) 

a/3 

where we have defined 

5 a/3 (pi) = $ a/3 (pi) - ($ a/3 (pi)) EE • (5.4) 

Our problem is to understand the double-sum term on the right-hand side of eq. (5.3). 

For a fixed value of pi, each |$ a ^(pi)| is given roughly by its RMS value in the 
eigenstate ensemble, which is A<& a/ g(pi) as given by eq. (3.13). We have already seen that 
A$ aa (pi) is extremely small, and A$ a/ j(pi) is not going to be any bigger when a / 
In fact, using the gaussian approximation of eq. (3.14) in eq. (3.13), we find 

[A$ a/3 ( Pl )] 2 ~ [A$ QQ ( Pl )] 2 exp[-m(C/ Q - U P ) 2 L 2 /Mn 2 U a ] . (5.5) 

We see that we will have A$ Q!j g(pi) <C A$ aQ (pi) if \U a — Up\/U a is much bigger than 
(h 2 /mUaL 2 ) 1 / 2 ~ A a /A^ 1 / 2 L, a very small number. [Note, though, that the precisely 
gaussian form of the falloff is an artifact of eq. (3.14).] For simplicity, let us assume that 
A/U < X/N X / 2 L, where A = (2nh 2 /mkT) 1 / 2 . Then for the range of a and (3 of interest, 
each A<E> a/ g(pi) is given by the right-hand-side of eq. (3.15) with U a replaced by U. We 
will need only the crudest approximations here, and so we write 

|$ a/3 ( Pl )| ~ A$ Q/3 ( Pl ) ~ {L/\y 3N/2 . (5.6) 

However, we expect that the phase of <& Q /?(pi) varies wildly with a and (5. 

Let A^c be the number of nonnegligible C a 's which appear in eq. (5.3); iVc can be 
defined precisely via iV^T 1 = ^ Q |C a | 4 , and crudely estimated as iVc ~ nA, where n 
is the energy level density near U = U [cf. eq. (4.6)]. The order of magnitude of each 
nonnegligible \C a \ is then N^ 1 ^ 2 (so that ^ Q |C a | 2 = 1). 

Now consider doing the double sum in eq. (5.3). If the phases of the C a 's are not 

carefully correlated with those of the <& Qi g(pi)'s, then each of the two sums will yield a 

l li 

random-walk result: the square root of the number of "steps", N c ' , times the size of each 
step, C a ~ N^ 1 ^ 2 . With an overall factor of (L/X) SN ^ 2 from eq. (5.6), we get 

J2c:C e^- u ^^( Pl ) ~ (L/\y 3N/2 , (5.7) 

which is again extremely small. If we let A/U be bigger than X/N 1 / 2 L : then the result is 
even smaller. 
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On the other hand, we can certainly set up an initial state which is very far from 
thermal. For example, we could give each particle the same initial value (to within quantum 
uncertainties) for its individual momentum, leading to an / QM (Pi)0) which is sharply 
peaked at that value. In this case, though, the phases of the C a 's must be correlated 
with those of the 5>a/?(pi)'s in exactly the right way to produce the desired nonthermal 
distribution / QM (pi,0). In this case, we want to see what happens as time evolves. 

Let us begin with eq. (5.3) at t = 0, with the phases of the C a 's carefully chosen to 
give us a nonthermal distribution. Now let the clock run. Each of the off-diagonal (a ^ f3) 
terms in the double sum begins acquiring an extra phase; there are roughly Nq off-diagonal 
terms in all. The growing phase of each individual term will cause its contribution to the 
sum to have a random orientation in the complex plane once \U a — Up\t/% > 2tt. We will 
say that such a term has decohered. The first terms to decohere (those with the largest 
difference between U a and Up) do so at a time t ~ h/A. The fraction of terms which 
have decohered at later times is given roughly by (A — h/t) 2 / A 2 . Thus the fraction of still 
coherent terms at this time is roughly 1 — (A — h/t) 2 / A 2 ~ (h/A) ft for t ^> h/A. 

Now, each of the coherent terms should give its usual contribution to the sum, since its 
phase is still (almost) properly aligned, but the net contribution of all the coherent terms 
will be suppressed by a factor of 0(h/ A)t _1 due to their reduced population. On the other 
hand, the terms which have decohered will contribute with random phases. Since almost 
all terms have decohered for t > ft/A, their total contribution will be given by eq. (5.6), 
and is negligibly small. Thus, overall we expect that any nonthermal features present in 
the initial distribution / QM (pi,0) will decay away with time like 0(h/ A)t -1 . 

We can test this conclusion with a very simple example. The system we will analyze 
is classically integrable, and so chaos plays no role in the following discussion. 

Consider a single particle with mass m = 100 in a two-dimensional circular box with 
radius R = 1; we also set h = 1. The initial wave function for the particle is taken to be 

V>(x,0) = 7r- 1 / 2 a- 1 exp(zp -x)exp(-x 2 /2a 2 ) . (5.8) 

This is a gaussian wave packet of width a at the center of the box, moving with momentum 
Po- If we Fourier transform into momentum space, we get 

^(p,0) = 7r- 1 / 2 aexp[-ia 2 (p- Po ) 2 ] . (5.9) 

We will take a = 0.1 and po = 100. Classically, the particle has energy E = p 2 /2m = 50, 
and just bounces back and forth, changing directions at times t = 1, 3, 5, ... . Quantum 
mechanically, the uncertainty in the particle's energy is A = po/^/2ma ~ 7. We can 
solve this problem exactly by expanding tp(p, 0) in the energy eigenstate basis and using 
eq. (3.4). Then we can compute the probability density for the particle to have its initial 
momentum at time t; that is, we compute |"i/>(po, t) | 2 . Let us see what we should expect 
for this quantity, based on the general arguments outlined above. 
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First of all, note that the infinite time average of jVKPcb t)\ 2 is given by 

lim -^eft |^( Po ,t)| 2 = X]|a| 2 |^(po)| 2 , (5.10) 

where V'a(Po) is an energy eigenfunction in momentum space. We expect | V'CPo? t) | 2 to 
approach its infinite time average at late times, with late-time fluctuations of the same 
order of magnitude. If | -0(po, t) | 2 does not start out close to its infinite time average at 
t = 0, then it should decay towards that value like 0(/i/A)t _1 . 

The exact result for |^(po, t) | 2 , normalized to its value at t = 0, is shown in fig. 1. We 
clearly see the classical bounces, as the probability to have the initial momentum drops to 
zero at the first turning point, t = 1, then returns at t = 3, etc. However, the quantum 
probability does not return to its initial value, but follows the (h/A)/t envelope predicted 
by the general argument outlined above. It finally drops down to its infinite time average, 
with fluctuations of the expected size. Thus we see that the simple phase decoherence 
argument works very well for this example. 

If we now compare the right-hand sides of eqs. (5.1) and eq. (5.10), we see the analogy 
between the infinite time average of [^(po, t) | 2 and the eigenstate ensemble average of 
/ QM (pi,t). The main difference is that, because we have integrated out almost all the 
degrees of freedom, the fluctuations of / QM (pi>*) about its average value are very small. 

There is more we can learn from our simple example, however. First, we have seen 
that the classical motion is reflected in the quantum probabilities, a fact which is expected 
to be true for classically chaotic systems as well [28]. This means that quantum initial 
states which can be understood as representing classical initial conditions may thermalize 
even faster, due to the effects of classical chaos. If we can follow classical trajectories (with 
initial quantum uncertainties) for some time, and classical chaos spreads these out over a 
constant energy surface in phase space, then the system has thermalized classically. This 
argument may be needed in cases where the initial distribution is so far from thermal that 
the quantum 0(/i/A)t _1 decay of its nonthermal features seems to take too long. It also 
shows that the 0(/i/A)t _1 rule need not be related to more traditional diffusion times, 
which are more likely to reflect classical physics. 

Second, consider fig. 2, which shows the |C a |'s for this problem plotted against the 
corresponding energy eigenvalues, E a . While they form a nice envelope, with mean energy 
E = 50 and uncertainty A = 7, there is a great deal of fine structure. This is needed to 
get the very special initial state, localized at the origin and moving in the x direction at 
a particular speed. We cannot, therefore, think of replacing C a by a smooth function of 
E a . This is unfortunate, since if we average eq. (4.4) over a smooth distribution of energy 
eigenvalues, the contributions of the long periodic orbits are suppressed [24,25,10], and the 
problem becomes much more tractable. However, as we see in fig. 2, such smoothing is 
physically far too restrictive, since it would prevent us from considering a wide variety of 
initial states which we could actually prepare in a real experiment. 
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Let us summarize the results of this section. If we construct a particular initial state 
for the hard-sphere gas at low density by superposing energy eigenstates, each of which 
individually satisfies Berry's conjecture, then we find that, at sufficiently late times, the 
quantum mechanical prediction for the probability that any one particle has momentum p 
is simply given by the Maxwell-Boltzmann distribution of eq. (1.1). The probability that 
this will not be the case is negligible, if we wait long enough. Any nonthermal features of 
the initial distribution for a single particle's momentum will decay away at least as fast 
as 0(h/ 'A)t _1 , where A is the quantum uncertainty in the total energy. Faster decays are 
possible, and likely if the initial state has a classical interpretation. Absolutely no averaging 
over initial states or times of measurement is needed, in contrast with the classical case. 

This concludes our analysis of the high energy, low density, hard-sphere gas in the case 
that we assume no symmetries of the wave function on exchange of individual particles. 
We will return to discuss lower energies and higher densities in Sect. VII. Now, though, we 
turn our attention to wave functions which are either completely symmetric or completely 
antisymmetric functions of the positions of the N particles. 

VI. BOSONS AND FERMIONS 

The detailed analysis in Sects. II and III required the assumptions of high energy 
(A a < a) and low density (Na 3 <C L 3 ), which combine to give A^A 3 <C L 3 . In Sect. II we 
noted that this is precisely the condition needed for the Bose-Einstein and Fermi-Dirac 
distributions to be well approximated by the Maxwell-Boltzmann distribution. Neverthe- 
less, even though the corrections due to quantum statistics may be numerically small, a 
valid formalism should be able to reproduce them. In this section we will see that the 
present formalism meets this test. 

Define a symmetrization operator V± via 

V±f( Pl ...p N ) = ^ ^(±l) p /( Pll ,...,p JJV ), (6.1) 

perms 

where the sum is over the permutations of the indices, and P is even (odd) if the permuta- 
tion is even (odd) . Now we can construct completely symmetric and antisymmetric energy 
eigenfunctions analogous to those of eq. (2.2): 

^(X) = N± J d 3N P A a (P) 5(P 2 - 2mU a ) V± exp(zP-X/ft) . (6.2) 

We will need the generalization of eq. (2.9). Using the same assumptions as in Sect. Ill, 
namely A Q < a and Na 3 <I 3 , we get 

CS*(P)^(P')) EE = <W (A^) 2 ^(P 2 - 2mU a ) V ± S 3 D N (P - P') . (6.3) 



17 



We now want to compute 

<^a(Pl)) EE = / ^ 2 ...^<^*(P)^(P)) ee . (6.4) 

For nonsymmetric wave functions, we found in Sect. Ill that ($ QQ (pi)) was equal to 
a Maxwell-Boltzmann distribution at a temperature T a related to the energy eigenvalue 
U a by the ideal gas formula U a = ^NkT a . We therefore expect to find that (^ , q Q (pi)) ei 
is given by a Bose-Einstein distribution f B E(Pi,T a ) = / + (pi,T Q ) or a Fermi-Dirac dis- 
tribution / FD (pi,T a ) = / _ (pi,T a ). In statistical mechanics, these are usually computed 
using the grand canonical ensemble, but in our case the number of particles is firmly fixed 
at N. Thus we expect to find that ($™(pi)) EE is equal to / ± (pi,T Q ) as given by the 
canonical, rather than the grand canonical, ensemble. Relevant formulae from the less 
familiar canonical ensemble are gathered in the Appendix. 

Let us warm up by computing N^- The normalization condition we need is 

j dV--.rf 3 Piv<^*(P)^(P)> EE = l. (6-5) 

We must first evaluate 

V± S 3 D N (P - P')| p , =p = ^ J2 ( ±:L ) P S *d(Pi - P»i) • • • S 3 d (pn ~Pi N ) , (6.6) 

perms 

where <5|>(p) is assumed to satisfy the substitution rules of eq. (2.8) with N — > 1. Now 
examine a particular term in this sum. If a particular momentum is paired with itself, we 
will say that it comprises a "1-cluster." If a particular momentum is not paired with itself, 
it will be part of an "m-cluster" of momenta which are all set equal to each other by the 
(approximate) delta functions. For each term, let m/ be the number of momenta in the 
Ith cluster, with m\ < . . . < mc- Also let C m be the number of m-clusters, and C be the 
total number of clusters. Obviously, we have the relations 

mi + ■ ■ ■ + mc = N , 

d + . . . + C N = C , (6.7) 
d + 2C 2 + . . . + NC N = N . 

Each term in the sum in eq. (6.6) can now be labeled by a set of nondecreasing integers 
{m} = {mi, . . . , mc}- The number of terms with the same label is 

c 

A {m} J](m,-l)!, (6.8) 
i=i 
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where (mi — 1)1 counts the number of ways momenta in the Zth cluster can be rearranged 
without breaking it into smaller clusters, and 

A{m} ~ (m 1 \...m c \)(C 1 \...C N \) (6 - 9) 

counts the number of inequivalent ways of assigning momenta to clusters. Each cluster has 
one redundant delta function, which results in a factor of (L/h) . Furthermore the Zth 
cluster contributes a factor of (±l) mi_1 to (±l) p . Now from eqs. (6.3) and (6.5) we see 
that to determine J\f^ we must multiply each term in eq. (6.6) by 5(P 2 - 2mU a ) and then 
integrate over all the momenta. Under the integral, when multiplied by a term labeled by 
{m}, we can make the replacement 

8{P 2 - 2mU a ) -> <y(m!P? + . . . + m c p 2 c - 2mU a ) , (6.10) 

All together, then, we have 

Wt) " 2 = h 3N £ A {m} {L/hf c I d* Pl ■ ■ ■ d 3 Pc 

{m} J 

c 

x 5(m lP j + ... + m c p 2 c - 2mU a ) ^(il)^- 1 ^ - 1)! . (6.11) 

i=i 

The sum is over all {m} with fixed N. We now make the change of variable p; — > m~ 1 ' 2 kj, 
which yields 

c 

[NtY 2 = h 3N J2A { m } (L/h) 3C I 3C (2mU a ) IJ(±l) m,_1 ("»i " l)!mr 3/2 , (6.12) 

{m} 1=1 

where Id(x) is defined in eq. (3.9). It turns out that terms with C ^> 1 dominate, and so 
we can use the large- C formula 

h 3N (L/h) 3C hc(2mU a ) ~ A^(L/A Q ) 3C h N (2mU a ) , (6.13) 

where A Q = (27r/z 2 /m/cT Q ,) 1 / 2 , to rewrite eq. (6.12) as 

c 

[Nt)~ 2 = \l N l3N(2mU a )J2A { m } II( L /Aa) 3 (il)™'" 1 ^, - 1)! mr 3/2 . (6.14) 

{m} 1=1 

Now we can apply the Mayer cluster-expansion theorem [29] , which can be written as 

c d N ( 



{m} 1=1 k L m=l 



(6.15) 

z=0 
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where, in our case, 

W m = (L/\ a f (±ir- 1 (m - l)!m" 3 / 2 . 
From eq. (A. 3) of the Appendix, we have 



m=l 



— W m = (L/Kfg± 2 (z). 



(6.16) 



(6.17) 



So putting all of this together, we find 



(A/;T 2 = X 3 a N I SN (2mU a ) ^\exp[(L/\ a y g f /2 (z) 



d 



N 



dz N 

= A 3JV I 3N (2mU a ) N\ Z± , 



i 3 „± 



z=0 



(6.18) 



where Zq is the canonical partition function for a gas of noninteracting bosons (+) or 
fermions (— ) at a temperature T a in a box of volume L 3 [cf. eq. (A. 8)]. 

Clearly we are on the right track! Now we have to do it all over again, this time 
leaving one of the N momenta unintegrated. 

Following the same logic which led us to eq. (6.12), we get 



c 



{m} 



i=l 



x (il)™*- 1 ^ -1)! JJCil)"*'- 1 ^ -l)!mf 3/2 



(6.19) 



The differences from eq. (6.12) arise as follows. First, we must choose which cluster contains 
the unintegrated momentum pi; this gives the sum over i = 1 to C. Then we must choose 
which of the rrij momenta in the ith cluster is unintegrated; this gives the factor of m^. 
Now we have overcounted by N, which results in the factor of 1/N. The change to the 
subscript and argument of I results from not integrating pi, and the factor of m~ 3 / 2 is 
missing because we did not have to rescale pi. 

Again, terms with C ^> 1 dominate, and so we have 



I 3C -3(2mU a - rriipj) ~ {2ixmkT a ) 3/2 exp(-mjp 2 /2mA;T a ) I 3C (2mU a ) . 
Then using eq. (6.13), we get 



(6.20) 



c 



fc(Pi)) EE = Wt) 2 X 3 f IsN(2mU a )J2A {m} J2 V ^U W ^ ' 

{m} i=l l^i 



(6.21) 
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where W m is given by eq. (6.16), and 

V k = (L 3 /Nh 3 ) (±l) k - 1 kl [exp(-p?/2mT Q )] fc . 
Starting with eq. (6.15), it is not hard to prove a generalization of it which reads 



(6.22) 



c 



E A i^} E V *n< n W ™i 

{m} i=l l=£i 



d N \ 

dz N { 



OO JU 

fc=i 



fc! 



exp 



E 

m=l 



ml 



(6.23) 



z=0 



In the present case, we have 

OO 



, , Nh 3 exp(p 2 1 /2mT a ) t z ' 

Combining eqs. (6.17), (6.18), (6.21), (6.23), and (6.24), we finally get 



(6.24) 



1 1 d 



N 



( $ aa(Pl)) EB = ^± w 



N\ dz N 



exp 



{L/xfgi /2 {z) 



Nh 3 exp(p^/2mT a ) =p 



= / ± (pi,T a ), 



z=0 

(6.25) 



where f^(pi,T a ) is the Bose-Einstein (+) or Fermi-Dirac (— ) distribution as predicted 
by the canonical ensemble [cf. eq. (A. 12)]. As expected, then, symmetrization or anti- 
symmetrization of the wave function changes the statistics from Maxwell-Boltzmann to 
Bose-Einstein or Fermi-Dirac . 

For the last time, we must study the fluctuations of <E> a a(Pi) that are predicted by the 
eigenstate ensemble. The relevant object is A$^(pi), defined by the obvious replacements 
in eq. (3.12). A$J /3 (p i ) is then given by eq. (3.13) with V± acting on S 3 D N (P-P')- Explicit 
evaluation of A$^(pi) is a fearsome combinatoric problem, but luckily a simple variation 
of the general argument presented after eq. (3.13) still applies, and can be used to show 
that A$^(pi) is very small compared to (3>a a (pi)) EB (^ > ^(Pi)) EE - Therefore eigenstate 
thermalization still holds, and the previous analysis (in Sect. V) of time evolution still 
applies. 



VII. DISCUSSION AND SPECULATION 



Let us begin with a brief recap of the central results. Berry's conjecture, as applied 
to a gas of iV hard spheres in a box, states that each energy eigenfunction appears to be 
a superposition of plane waves with wavelength fixed by the energy eigenvalue, but with 
random phases and gaussian random amplitudes. It is expected to apply only to systems 
which are classically chaotic, and has been found to be valid (with corrections that do 
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not affect our conclusions) in simple chaotic systems. Given Berry's conjecture for the 
hard-sphere gas, we have discovered the phenomenon of eigenstate thermalization: each 
energy eigenstate predicts a thermal distribution for the momentum of each constituent 
particle. This distribution is Maxwell-Boltzmann, Bose-Einstein, or Fermi-Dirac, de- 
pending on whether the energy eigenfunctions are nonsymmetric, completely symmetric, 
or completely antisymmetric functions of the N particle positions. Then, a superposi- 
tion of energy eigenstates with a small fractional uncertainty in the total energy will also 
appear to be thermal, unless the amplitudes and phases of the superposition coefficients 
are carefully selected to avoid thermal behavior. If this is done initially, then the usual 
phase changes produced by hamiltonian time evolution destroy the needed coherence, and 
any nonthermal features disappear as (9(/i/A)t -1 , where A is the uncertainty in the total 
energy. However, classical effects which are reflected in the quantum theory can result in 
faster thermalization. 

All of the analysis in Sects. II, III, and VI was done in the limits of low density: 
Na 3 <C L 3 , where a is the radius of a hard sphere and L s is the volume of the box, and high 
energy: A Q < a, where A a = (2nh 2 /mkTa) 1 / 2 is the typical wavelength of one particle when 
the energy eigenvalue is U a = ^NkT a ; numerically this means T a > (300 /ma 2 ) Kelvin, 
where a is in angstroms and m is in amu. But what happens if we relax these constraints? 

There are no fundamental difficulties with carrying out the analysis for moderately 
higher densities. All we need to do is use the exact formula for the smeared delta function in 
momentum space, eq. (2.7). In practice, though, this greatly complicates the calculations. 
It would be very interesting to try to develop some sort of perturbative (in a/N 1 ^L) 
analysis, and compare the results with more standard treatments of the hard-sphere Bose 
or Fermi gas [30]. 

Lower energies present an entirely different problem, since if we go low enough in 
energy, Berry's conjecture will break down. The question is, how low can we go? The 
generic expectation is that Berry's conjecture will be valid if the relevant wavelengths are 
small enough to "see" the features which produce classical chaos [16]. For the hard-sphere 
gas, the relevant feature is the nonzero radius a of each particle, which leads to A Q < a. 
However, this may not be good enough at high density [31]. Classically, if the density 
is large enough to result in very slow diffusion of the particles, then their positions will 
be correlated over long times; the Lyapunov exponents are all very small. We would then 
naturally expect that these correlations are reflected in the quantum energy eigenfunctions, 
which would mean that Berry's conjecture is not valid. In this case, a possible alternative 
criterion for the validity of Berry's conjecture is A Q < £, where £ is the classical mean free 
path of a particle, which can be much less than a. 

Whatever the correct criterion turns out to be, at a low enough energy Berry's con- 
jecture will break down, and we must ask what happens at lower energies. One possibility 
is that eigenstate thermalization will still be valid for a wide range of energy eigenvalues, 
even though Berry's conjecture is not. The reason for this speculation appears in the re- 
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suits of eqs. (3.15) and Sect. IV. In eq. (3.15), we see that the fluctuations about the mean, 
thermal value of $ aa (pi) in the eigenstate ensemble are extremely small; experimentally, 
we can tolerate much larger fluctuations. Thus, we may also be able to tolerate significant 
violations of Berry's conjecture without destroying eigenstate thermalization. In Sect. IV, 
this speculation receives some more support. Scars represent violations of Berry's conjec- 
ture which are quite obvious when one looks at the Wigner density of an energy eigenstate 
in phase space, since there the scars appear with a "signal-to-noise" ratio of 1:1. Once 
we integrate out all of the coordinates and most of the momenta, however, the scars fade 
away almost completely. The same should be true of more generic violations of Berry's 
conjecture. Thus, eigenstate thermalization may still be valid at energies well below the 
threshold for the validity of Berry's conjecture. 

If we go even lower in energy, though, presumably eigenstate thermalization will even- 
tually cease to be valid. If the system is this low in energy, it will not be able to thermalize 
itself. To find thermal behavior in a system below its threshold for eigenstate thermaliza- 
tion, we must couple it to an external heat bath, such as the refrigeration apparatus in a 
low-temperature experiment. Of course, once we have contact with a large, pre-existing 
heat bath, all the usual results of statistical mechanics can be applied without further 
worry. 

The basic question we have been trying to address in this paper is how such a heat 
bath might form in the first place. We have seen that this will happen for a hard-sphere gas, 
provided that Berry's conjecture is satisfied by the energy eigenstates which are superposed 
to form the initial state. Whether or not other mechanisms exist for self-thermalization of 
isolated quantum systems is an open question, one to which we hope to return. Meanwhile 
we believe that the present results constitute a new foundation for quantum statistical 
mechanics. In particular, we have at least one answer to the question of which quantum 
systems will approach thermal equilibrium. It is satisfying that this answer (those systems 
which obey Berry's conjecture) is closely related to the answer from classical physics (those 
systems which exhibit chaos). In fact the situation in the quantum theory is even better 
than it is in the classical theory, because we no longer need to consider an ensemble of initial 
states. Each and every superposition of energy eigenstates obeying Berry's conjecture 
will eventually yield a thermal distribution for the momentum of a constituent particle, 
provided that we wait long enough. Absolutely no averaging of any kind is needed: not 
over initial states, not over times of measurement, and not over hamiltonians. 

Finally we would like to comment on the much-discussed question of an appropriate 
definition for quantum chaos. Some time ago, van Kampen suggested that quantum chaos 
be defined as "that property that causes a quantum system to behave statistically" [32]. 
If we replace "behave statistically" with "obey the laws of statistical mechanics," then we 
have seen that the key feature is Berry's conjectured properties of the energy eigenstates. 
In particular, properties of the energy eigenvalues (such as GOE rather than Poisson 
statistics for the unfolded level spacings [33]) have played no role at all in the present 
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work. Steiner has suggested [21] that Berry's conjecture be elevated to the status of the 
best definition of quantum chaos, a proposal which we see to be equivalent to (our version 
of) van Kampen's. More generally, in quantum mechanics, where time evolution is always 
linear and therefore essentially trivial, the only place to encode the complexities of the 
classical limit is in the energy eigenfunctions: that is where quantum chaos, like thermal 
behavior, must be sought. 

I am grateful to Walter Kohn, Jim Langer, Joe Polchinski, Andy Strominger, and es- 
pecially Doug Scalapino for helpful discussions. I also thank Matthew Foulkes for pointing 
out an error in the original version of this paper, and the anonymous referee for several 
suggested improvements. This work was supported in part by NSF Grant PHY-91-16964. 



APPENDIX. THE CANONICAL ENSEMBLE 
FOR BOSONS AND FERMIONS 



We will use a notation close to that of [5]. We ignore spin degrees of freedom. 

The canonical partition function for N noninteracting bosons or fermions in a box is 
given by 

oo N± 

Z C = II E e ~ PUlEl <^,n 1+ n 2 + ... , (A.l) 
1=0 n,=0 

where (3 = 1/kT, E\ is the Zth energy eigenvalue for a single particle in the box, and 
N + = oo for bosons and N_ = 1 for fermions. 

Introducing the fugacity z, the grand canonical partition function for N noninteracting 
bosons or fermions in a box is given by 



N± 



Z GC = U 



ni -f3niEi 



1=0 m=0 



exp 



= exp 



= exp 



T $>g(l T ,e-^ 

1=0 

(va) 3 4 2 (. 



(A.2) 



In the third line, we have replaced the sum over levels by an integral over momenta 
(without separating out the zero mode, which would be necessary for a discussion of 
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Bose condensation). In the fourth line, we have introduced the thermal wavelength A 
{2nfi 2 /mkT) 1 / 2 , and the Lerch transcendents 



m — l^m 



(A.3) 



Z GC must be supplemented with the condition 

N = Z §- z l ° gZ GC 
= (L/\fg± 2 {z) 



(A.4) 



which can be thought of as fixing the value of z. We will call the positive real solution of 
this equation zq. 

The relation between Zq and Zq C can be found by writing the Kronecker delta in 
eq. (A.l) as 

(A.5) 



>N,n 1 +n 2 + ... — 



2ni 



dz z - N ~ 1 z ^+^+- 



where the contour encloses the origin. Substituting this into eq. (A.l) and using the first 
and fourth lines of eq. (A. 2) yields 



1 

2tH 

1 

2n~i 



dzz- N ~ x Z± 



GC 



dz z 



-N-l 



exp 



(L/X) 3 9^ /2 (z) 



(A.6) 



Evaluating this integral approximately by stationary phase, treating both iV and (L/X)' 
as large, results in 

.-1/2 ^ 

(A.7) 



2tt (L/\) 3 g± 2 (z ) 



where Zq C is to be evaluated at z = zq. The fractional error in this approximate equality 
is of order 1/N. Note also that, using Cauchy's theorem, we can rewrite eq. (A.6) as 



7 ± 

c N\ 8z N 



exp 



(A.8) 



z=0 



We would now like to compute the expected fraction (f ± (p))d 3 p of particles with 
momentum in a range d 3 p around p. The expected fraction of particles with energy E\ is 
given in either formalism by 



(ft) = (nf)/N 



1 1 dZ ± 



N(3 Z± dEi ' 
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(A.9) 



where Z^- is either or Zq C . Then converting to the normalization required for con- 
tinuous momenta gives (/ ± (p)) = (L//i) 3 (/ z ± ) with E\ = p 2 /2m. In the grand canonical 
case this gives the well-known result 



</ ± (p)) GC Nh S e P P y2m TZo 

In the canonical case, this procedure gives 



(A.10) 



c-^^-Jdz z~ N ^ exp 



{L/\fgi /2 {z) 



L 3 



(A.ll) 



Approximate evaluation of this integral by stationary phase gives {f^(p)} c = (/ ± (p))gc' 
again with a fractional error of order 1/N. Also, we can again use Cauchy's theorem to 
write 



1 1 f 

</ ± (p)>c = ^±m d^r |ex P [(L/A) 3 4 2 (,)j ^ ePp2/ ; mTz 



(A.12) 



z=0 



In the main text, we simplify the notation a bit via (/ ± (p)) c — > / ± (p, T). 
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FIGURE CAPTIONS 

Fig. 1. Solid line: |V>(p , t) | 2 /|^(Po, 0) | 2 vs t for a single particle in a two-dimensional circular 
box; the initial wave function is a narrow gaussian at the center with momentum po- 

Classically, the particle bounces off the wall at t = 1, 3, 5, Dashed line: (2ttJi/ A)/t, 

where A is the uncertainty in the energy. Dotted line: the infinite time average of the 
solid line. 

Fig. 2. Expanding the wave function of fig. 1 in energy eigenstates yields expansion coeffi- 
cients C a ; here \C a \ is plotted vs the energy eigenvalues E a . There are 1736 energy 
eigenvalues in the plotted range, 20 < E < 90. 
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